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ABSTRACT 

We analyse FORS2/VLT /-band imaging data to monitor the motions of both com¬ 
ponents in the nearest known binary brown dwarf WISE J104915.57-531906.lAB 
(LUH 16) over one year. The astrometry is dominated by parallax and proper mo¬ 
tion, but with a precision of ^0.2 milli-arcsecond per epoch we accurately measure 
the relative position change caused by the orbital motion of the pair. This allows us 
to directly measure a mass ratio of g = 0.78 ± 0.10 for this system. We also search for 
the signature of a planetary-mass companion around either of the A and B component 
and exclude at 3-cr the presence of planets with masses larger than 2 Mj^p and orbital 
periods of 20-300 d. We update the parallax of LUH 16 to 500.51 ± 0.11 mas, i.e. just 
within 2 pc. This study yields the first direct constraint on the mass ratio of LUH 16 
and shows that the system does not harbour any close-in giant planets. 

Key words: astrometry - brown dwarfs - binaries: visual - parallaxes - stars: indi¬ 
vidual: WISE J104915.57-531906.1 


1 INTRODUCTION 

The binary brown dwarf WISE J104915.57-531906.1 
(**LUH16), hereafter LUH 16 (Luhman 2013) represents a 
unique opportunity to study the properties of substellar ob¬ 
jects because of its proximity of 2 pc to the Sun. The system 
consists of a primary with spectral type L7.5 (component A) 
and a TO.5 secondary (component B) (Burgasser et al. 2013). 
Both show strong Li I absorption in their optical spectra 
(Faherty et al. 2014), which conhrms that they are the near¬ 
est brown dwarfs yet discovered. The system’s orbital period 
was estimated at 25-30 years (Luhman 2013; Mamajek 2013) 
and measuring the orbital motion offers the opportunity of 
constraining the component masses, hence determining yet 
unknown fundamental parameters of the system. LUH 16 is 
also an excellent target to search for orbiting planets (Gillon 
et al. 2013; Boffin et al. 2014), allowing us to investigate 
planet formation around brown dwarfs. 

2 DATA REDUCTION 

We analysed FORS2 images of LUH 16 taken between April 
2013 and May 2014 that were retrieved from the ESO archive 


(Programmes 291.C-5004 and 593.C-0314; PI: Boffin). The 
instrument setup and observation strategy used in these pro¬ 
grammes (/-Bessel hlter, several dithered frames per epoch, 
target position on CCD chipl, constraints on airmass and 
atmospheric conditions) are very similar to those of our exo¬ 
planet search survey (Sahlmann et al. 2014) and we reduced 
the data with our methods developed for that purpose (La- 
zorenko et al. 2009, 2014). Table 1 summarises the data we 
used and lists the epoch number, the mean date of the epoch 
exposures, the average airmass and the average FWHM mea¬ 
sured for star images. There are 22 epochs spanning 399 d, 
and every epoch consists of 16 to 42 usable individual expo¬ 
sures (Af) taken over At = 0.5 h on average, resulting in a 
total of 583 exposures. Figure 1 shows an example image. 

We used 585 reference stars located within a radius of 
2(2 of LUH 16 to measure the motion of the binary relative 
to the background held. They span a magnitude range of 
/ ~ 15 — 21, where 41 stars are within ±1 mag of LUH 16A. 
The components A and B were resolved in the images and 
we obtained their individual astrometry with an average per- 
epoch precision of 0.15 and 0.17 milli-arcsecond (mas), re¬ 
spectively. 


* Based on data obtained from the ESO Science Archive Facility 3 MASS RATIO MEASUREMENT 
under programme IDs 291.C-5004 and 593.C-0314. 

t E-mail; Johannes.Sahlmann@esa.int The standard model for FORS2 astrometry has seven free 

X ESA Research Fellow parameters (AaJ, Ma*, /i 5 , tu, p, and d) and reproduces 
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Figure 1. FORS2 image of LUH 16 in very good seeing conditions 
where the A and B components are well resolved. The image size 
is ~20" X 20", i.e. only a small fraction of the 4' X 4' field of 
view, with north up and east left. The star to the south-east is 
USNO-Bl.O 0366-0244607 {I = 15.7). 


Table 1. FORS2 data used in the astrometric analysis. 


No 

Mean date 
(UT) 

N{ 

At 

(h) 

Air- 

mass 

FWHM 

(") 

1 

2013-04-14T23:52:16 

16 

0.26 

1.25 

0.69 

2 

2013-04-20T01:28:07 

21 

0.31 

1.14 

0.66 

3 

2013-04-26T00:58:26 

21 

0.26 

1.14 

0.56 

4 

2013-05-06T00:46:31 

17 

0.36 

1.14 

0.83 

5 

2013-05-12T02;05:37 

21 

0.27 

1.23 

0.81 

6 

2013-05-20T23:26:08 

30 

0.76 

1.14 

0.78 

7 

2013-05-26T00:45:33 

19 

0.26 

1.20 

0.58 

8 

2013-05-31T23:48:51 

21 

0.26 

1.16 

0.68 

9 

2013-06-04T23:57:18 

21 

0.27 

1.19 

0.66 

10 

2013-06-09T23:18:32 

41 

0.58 

1.17 

0.69 

11 

2013-06-13T23:56:37 

34 

0.53 

1.24 

0.70 

12 

2013-06-16T23:23:53 

20 

0.32 

1.20 

0.71 

13 

2013-06-22T23:13:02 

24 

0.43 

1.22 

0.85 

14 

2014-02-05T06:00:02 

26 

0.67 

1.15 

0.74 

15 

2014-02-14T04:32:37 

42 

0.64 

1.19 

0.52 

16 

2014-03-10T03:31:31 

18 

0.61 

1.16 

0.68 

17 

2014-03-19T04;21:07 

34 

0.66 

1.15 

0.77 

18 

2014-03-30T02:57:14 

32 

0.67 

1.14 

0.75 

19 

2014-04-10T02:45:09 

32 

0.67 

1.15 

0.68 

20 

2014-04-25T23:58:25 

30 

0.67 

1.18 

0.66 

21 

2014-05-05T23:44:34 

29 

0.65 

1.16 

0.79 

22 

2014-05-18T23:35:28 

34 

0.67 

1.14 

0.75 


the astrometric measurements of a target ^ and Sm in RA 
and Dec, respectively, in frame m at time tm relative to the 
reference frame of background stars (Lazorenko et al. 2011; 
Sahlmann et al. 2014): 

CTm ~ AoiQ -f Pa* tm “t“ ri7 Ilajm P fl,x,m d f2,x,m \ 

dm = T PS im + -\-p /l,y,m + d /2,y,m, 

where AaJ,A(5o are coordinate offsets, pa*,ps are relative 
proper motions, and the parallactic motion is expressed as 
the product of relative parallax zo and the parallax factors 

^ We use the notation a* = a cos <5 throughout the text. 



Figure 2. Position change of LUH 16 B (open circles) relative to 
LUH 16A (filled red circle) over the course of the observations, 
corresponding to a section of the system’s orbit, which has to be 
highly inclined or very eccentric. Uncertainties are much smaller 
than the symbol sizes. 

Ila, . The parameters p and d model differential chromatic 

refraction (DCR) and their coefficients /i ,2 are functions of 
zenith angle, temperature, and pressure (Lazorenko 2006; 
Sahlmann et al. 2013). 

The epoch-averaged astrometry of LUH 16A and B is 
given in Table 1. Transformation to the ICRF was performed 
using 118 reference stars that are included in the USNO-B 
catalogue Monet et al. (2003). The residual dispersion of this 
mapping leads to an uncertainty of 80 mas in the absolute co¬ 
ordinates of LUH 16. However, all our results were obtained 
by fitting the 583 single-frame astrometric measurements, 
which are available upon request. 

3.1 Individual astrometry 

We first modelled the individual astrometry a\,SA and 
OBifc of LUH 16 A and B separately using Eq. (1). We aban¬ 
doned this approach because the resulting proper motions 
differed significantly between A and B by —212.1 ± 0.2 and 

322.1 ±0.2 mas/yr in pa* and pa, respectively, which is un¬ 
physical for a gravitationally bound system. In addition, the 
epoch residuals for A (0.9 mas r.m.s.) and B (1.1 mas r.m.s.) 
were anti-correlated and much larger than our measurement 
uncertainties, indicating that the motions of A and B are 
not independent and should be modelled globally. 

3.2 Orbital motion 

Figure 2 shows the position of component B relative to A 
obtained as a^ — a\ and 5 b — 5a for every observation epoch. 
Over the course of the observations, the relative separation 
decreases by 0"41 whereas the position angle decreased only 
slightly by 4.9° from its initial value of 314°. We thus detect 
the effect of the orbital motion of the system, with an am¬ 
plitude that agrees with the prediction of Burgasser et al. 
(2013) and that is consistent with the ~l"5 separation mea¬ 
surement of Luhman (2013). The direction of the relative 
motion indicates a highly inclined or very eccentric orbit. 
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3.3 Barycentre motion and mass ratio 

In a binary system, the astrometric orbital motions of the 
primary and the secondary relative to the barycentre are 
determined by the orbital period P, inclination relative to 
the sky plane i, eccentricity e, argument uj and time To of 
periastron, ascending node fl, and the component masses 
rriA and ms (e.g. Hilditch 2001). The latter are related to 
the barycentric semimajor axes aA and ds via dA/ds = 
ms/mA = q, where q is the mass ratio. Because the orbital 
parameters of A and B only differ by 180° in the argument 
of periastron, the position a!y,5j of the barycentre at any 
time is tied to the positions of the individual components: 

“ 7 = T^{<^A + qoiB) 

S,= jf-^iSA + qSB). 

We thus computed the barycentre positions of LUH 16 for 
a range of mass ratios according to Eq. (2) and fitted them 
with the model Eq. (1). The covariance matrices associated 
with the measurements of A and B were accounted for ac¬ 
cordingly and a linear least-squares fit (cf. Sahlmann et al. 
2014) was used to obtain the model parameters and their 
uncertainties for every value of q. 

Figure 3 shows the result for the mass ratio producing 
the smallest epoch The residual RMS is 0.23 mas, i.e. 
much smaller than for the individual fits discussed in Sect. 
3.1. However, it is still larger than the average epoch preci¬ 
sion of 0.11 mas. We attribute this to the decreasing relative 
separation of the pair, which in the second season approaches 
the median FWHM value of 0”7. As a consequence, the im¬ 
ages of A and B overlap more and more, which introduces 
errors of a few milli-pixel in our photocentre estimation. 

To determine the mass ratio of LUH 16, we inspected 
the epoch metric shown in Fig. 4, which has a parabolic 
shape with a well-defined minimum. The is minimum for 
q = O.TSlQ'Jg, where we derived the l-cr confidence interval 
from the g-values where the curve reaches a = 4-9.3039, 
appropriate for a fit with 8 free parameters. Table 2 shows 
the best-fit parameters at optimal and their systematic 
variation with q within its 1 -cr interval. 

In addition to the minimum x^ criterion, we also in¬ 
spected the correlations between the mass ratio value and 
the derived proper motions. Because proper motions have 
been measured with high relative precision from data span¬ 
ning 33 years (Luhman 2013), we can use them to determine 
which 5 -values are meaningful^. As shown in Fig. 5, there 
is good agreement between the q values determined via x^ 
and the range of values that are compatible with the proper 
motions of Luhman (2013). 

3.4 Updated parallax 

The analysis above also allows us to determine the parallax 
of LUH 16 with high precision. To convert the relative paral¬ 
lax in Table 2 to absolute, we used the methods described in 
Sahlmann et al. (2014, 2013) and the Galaxy model of Robin 
et al. (2003) to estimate a parallax correction of 0.28 ± 0.01 

^ We neglected the difference between absolute and relative 
proper motions, because we estimated that it is smaller than 1 
mas/yr for our FORS2 measurements. 




Figure 3. Top: The sky-projected barycentre motion of LUH 16 
measured with FORS2. Individual epoch measurements are shown 
with black circles and the model corresponding to the best-fit 
mass ratio of q = 0.78 is shown by the curve. The arrow indicates 
the proper motion per year. Bottom: Epoch residuals in RA (grey 
symbols) and Dec (black symbols) of the 7-parameter fit. 


Table 2. Astrometric parameters of LUH 16 determined solely 
by FORS2 data. The last two columns show the variation that 
corresponds to q = ± Oq. 


Param 

Unit 

Formal value 

Interval 

q 




0.78 

0.69 

0.88 

Aa* 

mas 

-1054.72 

± 

0.05 

±2.85 

-2.83 

ASo 

mas 

389.24 

± 

0.08 

±4.53 

-4.49 

tUrel 

mas 

500.23 

± 

0.06 

±0.11 

-0.11 

f^OL* 

mas/yr 

-2754.77 

± 

0.09 

±6.25 

-6.20 


mas/yr 

358.72 

± 

0.09 

±9.49 

-9.41 

PeS 

mas 

35.23 

± 

1.90 

±0.13 

-0.13 

C^efF 

mas 

-48.71 

± 

1.56 

±0.18 

-0.17 


mas, using the measured parallaxes of 175 reference stars. 
Taking into account the mass ratio uncertainty, we thus up¬ 
date the absolute parallax of LUH 16 to 500.51 ±0.11 mas, 
which agrees with the estimate of Luhman (2013) and cor¬ 
responds to a distance of 1.9980 ± 0.0004 pc. 

3.5 Prospect for measuring individual masses 

Whereas currently the estimation of individual masses re¬ 
lies on theoretical models, it will be possible to determine 
them accurately by mapping the relative orbit over about 
half of a revolution. The physical relative semimajor axis is 
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Mass ratio q 

Figure 4. Epoch (solid curve) and residual r.m.s. (dashed 
curve) in micro-arcsecond (/ras) of the barycentre astrometry fit 
as a function of mass ratio. The solid vertical line indicates the 
minimum value and the dashed lines indicate the 1-a confi¬ 
dence intervals determined by = 9.3. 



Mass ratio g Mass ratio g 

Figure 5. Correlations between q and proper motion in RA (left) 
and Dec (right). Blue symbols with uncertainties show the results 
of the linear fit. Horizontal lines indicate the values and uncer¬ 
tainties given by Luhman (2013) and vertical lines indicate the 
adopted q value. 

proportional to Ma + Mb, which in combination with our 
mass ratio measurement yields direct individual masses. We 
attempted to constrain the relative orbit using the small 
arc measured with FORS2, by considering only two extreme 
configurations: a face-on eccentric orbit and an inclined cir¬ 
cular orbit. A face-on orbit is highly unlikely because it re¬ 
quires extreme eccentricities (> 0.95) and very short orbital 
periods (< 8 years) to come close to reproduce the mea¬ 
surements. On the other hand, a circular orbit provides an 
acceptable fit to the data (residual epoch r.m.s. of ~1.6 mas) 
for an inclination of 95° and an orbital period of 45 years, 
which is slightly longer than previous period estimates (Ma- 
majek 2013; Burgasser et al. 2013). The configuration shown 
in Fig. 6 is the best fit to the data with 4 free parameters 
(Ma = 0.06 Mq, q = 0.78, e = 0, w = 180° were fixed) and is 
compatible with the elongation of the ESO optical image in 
1984 discussed by Mamajek (2013). A highly inclined orbit 
is the most likely to be observed for a ensemble of randomly 
oriented orbits, however, this simplified and tentative orbit 
will have to be revised with future measurements. If the 45 
year period is accurate, we estimate that reliable individual 
masses can be obtained around the year 2030. 



Offset in Right Ascension (mas) 

Figure 6. Possible relative orbit of LUH 16 with a period of 45 
y, derived under the assumption that is is circular. 


4 LIMITS ON A TERTIARY COMPONENT 

The small residuals of the barycentre fit presented in the 
previous section assert that our model represents the data 
well. We can thus use the data to derive limits on a possible 
third body in the system, which would introduce an addi¬ 
tional astrometric signal. We first searched for significant 
signals in the residual periodogram, finding none. At a fixed 
distance, the astrometric signature of an unseen companion 
depends on its mass, the orbital period, and the geomet¬ 
ric configuration of the system. Accounting for the actual 
FORS2 sampling and the accuracy of our measurements, 
we then derived companion detection limits as a function of 
period and mass. 

We followed the procedure of Sahlmann et al. (2014), 
which consists in simulating 5000 circular orbits for each 
point in a grid of 100 periods (20-1000 d) and 60 compan¬ 
ion masses (0.1-5 Mj), i.e. 30 million pseudo-orbits, with 
randomly assigned parameters fl, Tq, i, the latter with a 
probability density function of |sini. For every orbit, the 
astrometric signal is computed at the sampling times and 
added to the observed astrometry, followed by the linear fit 
using Eq. 1. A companion of given mass and period is con¬ 
sidered detected if 99.7% of the corresponding pseudo-orbits 
have a residual r.m.s. larger than the one observed for the 
target. This procedure yields 3-cr mass detection limits as a 
function of period. 

We run these simulations for both components of 
LUH 16, where an additional step is required because the 
reference residual r.m.s. was obtained for the barycentre. A 
photocentre shift of component A and B will result in an ap¬ 
parent barycentre shift detected by our measurements that is 
smaller by a factor of l/(l-|-g) and q/{l-\-q), respectively (Eq. 
2 ). Therefore we consider host masses of Ma = 0.06 Mq, ap¬ 
propriate for an age of ~3 Gyr (Faherty et al. 2014), and 
Mb = qMa — 0.047 Mq and increased the reference r.m.s. 
by the respective reciprocal factor, which deteriorates the 
detection limit. We repeated this with Ma = 0.03 Mq, cor¬ 
responding to a younger age of ~0.5 Gyr. 

Figure 7 shows the maximum mass of any putative com¬ 
panion to LUH 16 A and B as a function of its orbital pe¬ 
riod. The sensitivity deteriorates for very short periods be¬ 
cause the astrometric signal decreases, at long periods be- 
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Figure 7. Companion mass limits for LUH 16 A (solid) and B 
(dashed) as a function of orbital period. Thick and thin lines 
correspond to primary masses of Ma = 0.03 Mq and 0.06 Mq, 
respectively. The maximum companion mass compatible with the 
measurements is shown, i.e. companions on and above the curve 
are excluded by the data at 3-cr confidence. 

cause those exceed the measurement timespan, and at 365 
d because of correlation with the parallax signal. However, 
for orbital periods between ~20 and ~300 days (~0.05-0.4 
au), the FORS2 data exclude at 3 -(j confidence the pres¬ 
ence of planetary companions with masses larger than ~2 
Jupiter masses (Mj) around A and B, assuming a primary 
mass corresponding a 3 Gyr age. At 100 d period, any com¬ 
panion more massive than Jupiter is excluded. If LUH 16 
is younger ~0.5 Gyr, the detection limits are a factor of 
two better. These limits also exclude the presence of a third 
brown dwarf in the system within the explored period range, 
except for a potential companion that has equal mass and 
7-band luminosity as the object it orbits, because for these 
systems the induced photocentre shift would be null. 

We thus excluded that either of the LUH 16 components 
harbours a giant planet more massive than twice the mass of 
Jupiter in a short period (20-300 d) orbit, which disagrees 
with the conclusions of Boffin et al. (2014). So far, the system 
appears to consist of only two components. 


5 CONCLUSIONS 

We analysed multi-epoch imaging observations of LUH 16 
and measured the astrometric motions of both the L- and T- 
dwarf components that show considerable orbital motion. By 
modelling the system’s barycentre, which follows the stan¬ 
dard parallax and proper motion, we determined the mass 
ratio of the two brown dwarfs to be g = O.TSlQ'Jg, which is a 
model-free measurement derived directly from the astrome¬ 
try and that does neither depend on the parallax value nor 
on the orbital parameters. In comparison with other ultra¬ 
cool binaries, this mass ratio is low but not exceptional (Liu 
et al. 2010). We also updated the parallax to 500.51 ± 0.11 
mas. 

As one of the closest stellar objects to the Sun, LUH 16 
is an ideal target for astrometric planet search. We used the 
FORS2 astrometry to search for a planet in the LUH 16 sys¬ 
tem and found no significant signal. We excluded the pres¬ 
ence of Jupiter-mass planets in short period orbits. Contin¬ 
ued monitoring with FORS2 or other instruments will im¬ 


prove the sensitivity to lower-mass planets and complemen¬ 
tary imaging, spectroscopic, and photometric observations 
will allow us to further investigate the presence of extrasolar 
planets in the immediate solar neighbourhood. The results 
presented here contribute to the detailed characterisation of 
LUH 16, which is a fundamental system for understanding 
brown dwarfs and their physics. Our mass ratio determi¬ 
nation helps with the interpretation of spectroscopic and 
photometric observations, leading to a better picture of the 
properties of LUH 16. 
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Table 1. Epoch astrometry of LUH 16A and B in the ICRF, after the effects of DCR have been removed. The quoted uncertainties 
correspond to the photocentre precision. The conversion to ICRF introduces an additional uncertainty of 80 mas in both RA and Dec, 
which was not incorporated here. 


Epoch 

(MJD) 

RA 

(deg) 

Primary 

ctra Dec 

(mas) (<i6g) 

<^Dec 

(mas) 

RA 

(deg) 

Secondary 
o-RA Dec 

(mas) (deg) 

<^Dec 

(mas) 

56396.9941 

162.31119262 

0.18 

-53.31833644 

0.23 

162.31071214 

0.23 

-53.31805513 

0.25 

56402.0614 

162.31116077 

0.16 

-53.31832531 

0.18 

162.31068269 

0.18 

-53.31804499 

0.18 

56408.0411 

162.31112537 

0.18 

-53.31831183 

0.20 

162.31064843 

0.21 

-53.31803272 

0.21 

56418.0342 

162.31107044 

0.19 

-53.31828831 

0.21 

162.31059589 

0.21 

-53.31801145 

0.23 

56424.0872 

162.31103953 

0.22 

-53.31827345 

0.25 

162.31056730 

0.25 

-53.31799796 

0.26 

56432.9754 

162.31099898 

0.18 

-53.31825165 

0.20 

162.31052781 

0.20 

-53.31797790 

0.21 

56438.0317 

162.31097762 

0.21 

-53.31823918 

0.24 

162.31050839 

0.23 

-53.31796661 

0.24 

56443.9928 

162.31095466 

0.22 

-53.31822436 

0.25 

162.31048650 

0.26 

-53.31795325 

0.26 

56447.9940 

162.31094034 

0.18 

-53.31821458 

0.20 

162.31047347 

0.20 

-53.31794448 

0.21 

56452.9718 

162.31092407 

0.15 

-53.31820274 

0.16 

162.31045826 

0.16 

-53.31793372 

0.17 

56456.9943 

162.31091192 

0.16 

-53.31819350 

0.18 

162.31044753 

0.18 

-53.31792530 

0.19 

56459.9747 

162.31090367 

0.17 

-53.31818664 

0.19 

162.31043962 

0.19 
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